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t^ . Abstract 
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We prove boundedness of gradients of solutions to quasilinear parabolic system, the main part of 

which is a generalization to p-Laplacian and its right hand side's growth depending on gradient 
^S| ' is not slower (and generally strictly faster) than p-l. This result may be seen as a generalization 

CN I to the classical notion of a controllable growth of right hand side, introduced by Campanato, 

over gradients of p-Laplacian-like systems. Energy estimates and nonlinear iteration procedure 
Q, I of a Moser type are cornerstones of the used method. 
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1. Introduction 

1.1. General statement of the problem. 

We are interested in obtaining a local boundedness of gradients of solutions to a following 



00 

(<— ^ , parabolic system in Q c 

(N 



■(A;,(Vm)),„, =/'(x,f,VM) i^\,..,N 



where the main part is a generalization of p-Laplacian and the right-hand-side grows as 1 + |Vm|" 
or IVmT with w specified further We say that a right-hand-side is a fast-growing one, when 

^ , w > p -I holds. 

C^_' The existing literature on the regularity issue of parabolic equations and systems is impressive. 

Let us recall that for equations the existing results are quite strong: even for the right-hand-side 
growth of 1 + Nu\'' one obtains C''" regularity of solutions: see classic monograph Ladyzhen- 
skaya et al. 11911 for the case p - 2 and DiBenedetto 181] for p e (1, oo). Many further general- 
izations are possible: for instance in Bartier and Souplet ^ the right-hand-side takes the form 
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e"|V«|'', which suffices for a boundedness of Vm. Moreover, this growth condition seems to be 
optimal, because there are blowup results for gradients of solutions to equations, the right-hand- 
sides of which grow faster than p - compare Souplet |24J]. In case of systems, the regularity 
results are much weaker One can construct irregular (i.e. unbounded or discontinous) functions, 
which solve homogenous paraboUc systems. For « > 2 it suffices for irregularity that the coeffi- 
cients A(x, t) of the main part are discontinous (and still bounded) or that there is a relevant non- 
diagonality of the main part - for details, consult Arkhipova [ij]. Nevertheless, there are many 
classes of main parts, which allow for higher regularity (even C'") in the homogenous case; 
these are: having structure close to Laplacian or p-Laplacian, like those studied in Ladyzhen- 
skaya et al. |l9ll or DiBenedetto [SO, respecti vely , or having main part depending solely on Vm: 



see well-known paper by Necas and Sverak Il22[l or more extensive research done by Choe and 



Bae ||7[]. As these papers consider homogenous systems, one may ask a natural question: what 
inhomogenous counterparts of such systems remain, in a certain sense, regular? The general 
answer is unknown, but there are several hints: on one hand, for the right-hand-side growing Uke 
1 -I- |Vm|''"' the regularity of the homogenous case seems to be retained - see DiBenedetto (SO; 
on the other, unlike for equations, one cannot have the right-hand-side growing as fast as |Vm|'' 
without further assumptions, even in case of a system with the simplest main part, i.e. an inho- 
mogenous heat system. Recall the classical counterexample: for « > 3 bounded but discontinous 
function ii(x) - A with unbounded weak derivatives solves u\ — t^ii' — m'|Vm|^ = (n - 1)t^ , 



for details - see Struwe 112511 . It turns out that in case of an inhomogenous system for p - 2 
one has to additionally assume a certain smallness in order to obtain regularity - for details, refer 
to Tolksdorf |!26y, Pingen pSy, Idone 11511 or even the classical Ladyzhenskaya et al. lll9ll . The 



regularity issue for a general nonlinear inhomogenous parabolic system with the right-hand-side 
growing at the rate 1 + IVmI"" for w possibly close to p, homogenous counterparts of which enjoy 
regularity, is not fully researched, especially for the case p ^ 2. There are several approaches 
to answering this question: some authors relax the notion of regularity by resorting to partial 
regularity - see for example classical papers of Italian school: Campanato l^, Giaq uinta and 
Struwe il3ll and newer ones: FanciuUo llllll . Frehse and Specovius-Neugebauer 11211 . 



^ _ Misawa 

I20I1 . Duzaar and Mingione JlOll : or by demanding a high integrability-type regularitjQ like in 
Naumann 12111 . Kinnunen and Lewis 01811 or Bensoussan and Frehse ||31] (in the last paper the 



growth of the right-hand-side may be polynomially arbitrarily large!). Certain systems with pe- 
culiar structure or two-dimensional ones (or at least close to them in some sense) enjoy also high 
regularity, even if they are much more general than a Stokes-type system; for results in this di- 
rection compare papers of Seregin, Arkhipova, Frehse, Kaplicky (and many others), for instance: 
Arkhipova |[l[], J. Naumann and Wolff lll6ll . Kaplicky 11711 . Zaj^czkowski and Seregin 12711 . 



In this note we focus on deriving a full regularity result, more precisely: the local boundedness 
of gradients, for a class of quasilinliear parabolic inhomogenous systems. Our goal is twofold: 
firstly to obtain results for a general inhomogenous parabolic system, the main part of which is 
analogous to the system considered in Choe and Bae ||7|], while retaining possibly general growth 
conditions for the right-hand-side. Secondly, to sharpen these results with respect to growth of 
the right-hand-side, restricting ourselves to less general systems, being close to /7-Laplacian. 



For similar result on the level of solutions, compare Giorgi and O'Leary 1114 1. 

Let us emphasize that we proceed in a manner typical for the regularity approach: we assume 

existence of solution m in a given class, which is often a deep problem itself, from which we 



^such results are especially interesting, as our result may be easily strenghtened via higher regularity 
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derive higher regularity. Moreover, we concentrate on a priori estimates while conducting the 
proofs; the rigorous version of computations is commented on in the conclusion. 

1.2. General definitions and assumptions. 
Consider parabolic problem in Q c K" 

u', - (A;,(Vm)).,„ = f{x, t,Vu) i^l,..,N (1) 

As all our results have a local character, any further specification of Q is irrelevant. 

We say that a vector valued function u 6 L"'(0, T; L~(D.)) n LP(0, T; W''''(f^)) is a weak solution 

todBiff 

I -u'd + AL(Wu)(p[ dxdt = I f'(x, t, Wii)(/)'dxdt V 

Globally, the following notions will be used; 

• (5g denotes the Kronecker delta, 

• Qr{xo, to) denotes paraboUc cylinder, i.e. Br{xo) x {to - R'\ to); when possible, cut short to 
Qr, 

• TjpR e C^{Qr) denotes standard parabolic cutoff function for Qp c Q^, when possible, 
cut short to T], which satisfies; 77 = 1 in Qp, t] - Q outside Qr, IV77I < c{R - p)^^, |;;_,| < 
c(R-p)-'\ 

Throughout the article, summation over repeated indices is in use. 

1.3. The structure of results. 

We show our results in the following order; 

1. First, we derive a general result for inhomogenous version of system analyzed in Choe 
and Bae \M, where ellipticity assumptions for the main part are generalized by introducing 
exponent q (Theorem |2]i. Here, loosely speaking, admissible growth for the right-hand- 
side is 1 + \Wu\''^\ so this result may be seen as parallel to DiBenedetto IS]]. 

2. Next, we admit faster growths for the right-hand-side, at the cost of assuming that the main 
part is closer to p-Laplacian, in the sense that it is not eruiched with terms involving q > p 
(Theorem|3]l. 

3. Finally, we state the result for the least general case, i.e. for the 3D p-Laplacian with the 
right-hand-side growing as |Vm|" (Theorem[l]i- 

Since the last result seems to be the most traceable one, let us give incentive to studying the 

technical remainder of this paper by stating Theorem [T]now; 

Consider u = («') E L~(0, T; L^{Q.)) n L'^(0, T; W''''(Q)) solving the 3D p-Laplace system 

u', - div(| VMr^VM') = fit, X, u, Vm) / = 1, 2, 3 

Theorem 1. Let Q be an arbitrary domain. Assume a growth condition: \f'(x, t, Vm)| < IVm]**^, w < 
p and initial integrabilityQ |Vm|^p < 00. If one of the following conditions is fullfilled: 



'From existence one has p = p, so this assumption may be void. It only helps to quantify the results when we have 
some additional knowledge on integrability. 
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1. w < p - 1, p = p; 

3. w < ^ and p > 2 — |p 
then Vu is bounded: 



IVmIz." < C|Vm| 






For the proof, see the end of the next section. 

Observe that for w > p-1 in the degenerate case (i.e. /? > 2) we have w > ^, so point 2. applies. 
In such case merely from existence, i.e for p - p, one has w < /? - |. This corresponds for p -2 
with Campanato's notion of controllable growth, stating that to obtain regularity, the growth of 
right-hand-side must be smaller than | (see for example Campanato [4]). Therefore, this result 
may be seen as a generalization of classical results over gradients of p-Laplacian-like systems. 
Utilizing results on high integrability of certain systems, one may relax the growth condition 
further For example, for the system analysed in Naumann [21] we have w < e + |, because 
Vm e U^'^{Q.j), which can be taken as the initial integrability L'' . 

In all our theorems there is no explicit assumption that w < p. In fact, the inequality w < p 
is enforced by the rigorous treatment of energy estimates. Simultaneously we know from the 
counterexample recalled in the introduction, that w = /;> is generally not admissible. Therefore, 
our results can be viewed as a way to quantify the possible boundedness of gradients by means 
of higher integrability. As in the just mentioned case of Theorem!!] for;? = 2 one has w < p-^, 
which can be boosted in some cases to w < e + |, because Vm e U^'^. For w - p one would need 
Vm e L^+5. 



2. Boundedness of gradient of solution 

As outlined in the introduction, first we prove the general theorem. As the cornerstone of the 
analysis is the energy method, we derive formal estimates for the sake of transparency. For 
a rigorous justification of the formal estimates please consult the conclusion of this note. We 
analyze solutions of 

mJ - (A;,(Vm)),., = f'{x, t,Vu) i^l,..,N (2) 

where the main part comes from Choe and Bae |7i] and the right-hand-side grows as 1 +|Vm|" , w < 
6 for a certain 6 < p, obtaining the boundedness of Vm. More precisely, one has the following; 

Theorem 2. Under the following assumptions: 
(AO) ellipticity-type: A[-^ is given by potential F e C^(M), F'{Q) > 0, as follows 

KiQ) = (F{\Q\)& 

and F enjoys ellipticity 

where: 

l<p<q<p + l<oo 
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(Al) growth-type 

where 
(A2) initial integrability 



\f'\<+c,\vur+c2 

p>w>Q, c; 6 L"'(Qr) 



Vm e m'^^iD.r), M := max(2, p,2q- p,w +\,2w - p + 2) 



with sq satisfying 



and 



np Mn 

sq>0, S0+2 + -^ >0 

2 2 



sq > p -2 for C2 i^O 

So > p - 2w - 2 for C2 -0 



gradient of solution to ^ is locally bounded; moreover, for any Qr^ C Q. with Rq < 1 the 
following inequality holds 



\Vu\ 



i."(G«o) 



C \ \^u\'°^^dxdt 



..o+2+T-^ 



+ C 



Proof. First we derive formal energy inequalities, then we implement an iteration scheme. 
Differentiate formally system (|2) to obtain 

u' - {A' , (Vm)m4,^).,, = (/'(x, t, Vm)),^ / = 1, .., a? 

' '11 

testing (|3]i by u'^ |Vm|'?7^ one gets 



(3) 



1 d 



s + 2dt Jb, 



\Vur-T]^dx 






- f fix, t, V«)[<.^,JV«r772 ^ ,„y v«|-2«*. ^.^„^^^;;2 + 2uijVu\^rjT^,^]dx 

•J Br 

+ ^ r |V«|''+2^77,«'x - r A' (V«y «* |V«|>,.,«fx (4) 
Consider /. Utilizing ellipticity assumption (AO) with ^' '■= u[ ^ one estimates the first summand 



of / as follows 



XfXp 



•J B[i fi 'J Bfi 



and because A[^ is given by potential F, from differentiation we estimate the second summand of 
/ 



I A' ,(Vm)m{„, 5|Vm|' ^u[ u'l , u'l i?dx ■ 



^"(|V«|)^!^+F'(|Vm|) 



^^'.A <A^ 



|Vm| |VmP 



«it,^|VMr^Mi^<,,M*,/7^£/x 



= . I /7V"(|V«|)|V«|-4(4„;/„i )(4„^„^ )^^+ 



•JBk ^ ' ^ ^ ' 

= \ (I V/(P),„ (I VhP),„ |V„P = 1 < (I V«p),„ /,i^ (I v«px,^ 

From the ellipticity assumption (AO) one has F"{\s\) > 0, F'{Qi) > therefore it holds: /^'(l^?]) > 0. 
This, in conjunction with the following computation: u'^ (|VmP).v„m\. (|Vm|^)j: < |VmP|V|VmPP - 
|VmP(|VmP)x^(|Vm|^)j;^, implies that equation ^ takes the form 



f A' iC^uA^ slVul'-^u' u\ ^ u\ rj^dx > 



(7) 



Summing up (|5]l and (|7]) we conclude that / satisfies 



>a\ IVmI'^^^IV^mI^IVmIV^jc (8) 
Inputting inequality (O into (JUi we arrive at 

f \VuY^^rf-dx + A \ IVm^^IV^mI^IVmI V-^^ < 

< c I [1 + |VMr''][(l + i)|V2M||VM|V + 2|VM|"+'?7|V77|]c/jt;+ 

where the last inequality is valid in view of growth (Al) and ellipticity (AO) assumptions. 
Absorb |V^m| from the right-hand-side of (|9]l using Young's inequality and integrate with respect 



to tune 

1 

s + 2 , Jb^ Jq^ 



sup ( \Vur^jfdx + iA-s) f \'7u\P-^\V^u\^\Vu\'T]^dxdt 



c r (l + |VM|'0|VM|'+'77|V77|£/xc/f+-^ f \Vu\'^^r]\ri,\dxdt+ 
Jq, s + 2 Jq^ 

c f \Vu\'[\Vu\'' + \Vu\^''-P]j]^\VT]\^dxdt + c{l + s) f r]-\Vun\Vu\^'''''*^ + \'^u\--P]dxdt (10) 

Jqr Jqr 

By estimates for derivatives of the cutoff function ;; we obtain 

^— sup r l^^urWdx + iA-s) f \Vu\''-^\V^u\^\'^u\'j]^dxdt < 
■5 + 2 ( Jbr jQg 

■ \ |Vm|'[|Vm|'' + |Vm| + |Vm|"+' + 



(2 + s)-^\Wuf + \Vuf + |Vm|2«"'' + (1 + s){\Vu\-''-''^^ + \'^u\--P)]dxdt (11) 

for Q < p < R < 1. Since for some w, p the exponents 2w - p + 2,2- p may be nonpositive, we 
estimate respective powers of |Vm| using |Vm|' as follows 

r |Vm|^[(1 + s){\Vu\^''-P^^ + \VutP)]dxdt <{l + s) f [1 + |VM|™''<'+2-p,.+2v.+2-rt]^^^^ 

Jqr Jqr 

(12) 

For the last inequality to hold, we must assume 

s > ma\{p - 2w - 2, p - 2) (13) 

Because summand |VmP"'' occurs only if C2 ^ in the growth condition (Al): |/'| < ci |Vm|*'' + C2, 
the above assumption ( fTlT l can be written as 

s > p — 2 for C2 ^ 

s > p — 2w — 2 for C2 — 

In the forthcoming iteration scheme we construct a growing sequence of i,, therefore it is suffi- 
cient to assume 

{So > p -2 for C2+Q 
So > p -2w -2 for C2 - 

which coincides with our initial integrability assumption (A2). 
By computation, the following inequality holds 

r |V(|VM|^77)|2(ixt/f < I IVuf^'l'^jjfdxdt + {p + sf f IVuf^'-^lV^ul^jfdxdt (14) 

Jqr Jqr Jqr 

Adding to both sides jjffvT Jq \^u\'''^''\^ri\-dxdt and considering properties of rj, as (s + p)~ < 
c(l + s^), we arrive from (ITTt . by virtue of (IT4l . at 

sup I |Vm|"V'^^+ r W{Wu\'^r])\-dxdt < "^ ^ j I 1 + IVMl'+'^t/xt/f (15) 

' J Br Jqr (R - P) JQr 
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taking into account (fT2l i if neccessary. Recall that by definition M - max(2, p,2q-p,w+l,2w- 

p + 2). 

By Holder and critical-Sobolev inequalities (respectively), one gets 






dt 



sup I |VM|'+^?7^iix 



' *J Bfi J *J t{)—R~ \_kJ Bn 

+ s') r 



(if < 



|V(|Vm|— 77)|^(ix^f < 



£(1 



1 + \Vu\'^'^dxdt 



(16) 



Inequality ( fT6b is our desired energy estimate, which we now iterate. Define recursively numbers 
sr- Si+i + M - p + Si + (si + 2)^, then 

n(p - M) 



n 



sq + n + 2 



n(p - M) 



Utilizing the initial-integrability assumption, i.e. io + 2 + « - ^^^ — - > 0, we have 

n(p - M) 



Si > oo: 



Si /^OO 



S(j +2 + n 



Let 



(1 + -„)' 

^i = i \Vu\"^^dxdt 
Jsg. 

then (fTSI l with rjR.^^ji. can be written as 



(17) 



|5r,+i|iA;+i 



C(l + 5f)(^| |5fi,|(l+.A/) 



/?:';f'A,.i < 



^P 



/2'+2\ 

c(i + /)i — 1 /?;.';2(i + 'A/) ^ «A,+i<c(i + /)^2'(i + ,A/)f (18) 

with p :- 1 + ^; the last inequality given by Ri :- -j^(l + 2^'). As we know from ( [TtI i that 
asymptotically i, behaves like/?', finally ( fTSl l folds to 

IA/+1 < C'lAf + C' 
which, by a standard computation (see Choe and Bae ||7|] for details), gives 

iA,+i<C''"'<"+(/ + l)C^'" 
From the above considerations one gets, using the definition of i/r 



(19) 






Qlir 



+" 1 ., (19} 

T'i+ — 



(C^"'^A^")^^+((/+l)C^'") 



i^oo so+2+n-^ 



^ Ci^r ~ + C (20) 



in view of (fTTI i. 

As Sj + M > oo, ( |20t in tandem with the initial integrability assumption gives the following 

uniform bound 



|VM|L»(e«„)<C f ivur^^'dxdt 



-Mr 
p 

+ c 



D 



In the next theorem we resign from the term possessing q > p in the ellipticity assumption. This 
allows us, in turn, to obtain bigger growths of the right-hand-side, as now it is possible to derive 
estimates for negative s > -j^. 

Theorem 3. Gradient of solution to ^ is locally bounded, under the following assumptions: 
(AO) ellipticity-type: A'^ is given by potential F € C^(M) as follows 

aUQ) = (^(lei)a 

and F enjoys ellipticity 

(Al) growth-type 

l/'l < c|Vm|"', p>w>0, ce L'°(Q.t) 

(A2) initial integrability 

Vm e L'^^^iilr), M ;= max(2, p,w + \,2w - p + 2) 

with sq satisfying 

So > max(--^, p -2w -2) 

moreover, for any Qr^ C Q with Ro < 1 following inequality holds 






C\f \Vu\'"^'^dxdt 

ve«„ ) 



+ c 



Proof. For s > Theorem |3] is a special case of Theorem |2] therefore it suffices to show it in 
the case of negative s. The only difference in the energy estimates is the lack of positivity of the 
left-hand-side term of ^, where sign of s plays a role 



r A' ,iyu)u{„, s\Vu\'-^u\ u'i ^ ui rj^dx 



it can be however estimated as follows 



s r TJ'F\\Vu\)\Vur%y,y_,^,u[u%ul,/x > sA f ri^Wur'-^W^fdx 

J Br *JBr 
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which allows for a following counterparty of (|8) 

JBg V 

From this inequality on, one proceeds identically as in proof of Theorem|2] D 

As a consequence of Theorem|3]we obtain Theorem[T] 

Proof of Tlieorem\J] The point 1. stems from theory in DiBenedetto ||8|] and it implies that in 
points 2. and 3. one can assume w > p - I- The rest of the proof follows from Theorem |3] 
Indeed, for p-Laplace system /I = A, so the ellipticity condition is fulfilled and the rest of the 
assumptions of Theorem|3]can be rewritten with so = p - M as follows: 

M:=max(2,2w-/7 + 2) 

p>maxlM- 1,^-^-2] (21) 

In the case w > j. M - 2w - p + 2 and ( 1211 1 gives w < ''^'^ for w < p - | and w < ''^ ^^ for 
w > p - |; as in view ofw>p-l and p > p the first condition is void, we obtain assumption 

^ p+4p-3 
W < ■ — J — . 

In the case w < j. M -2 and (1211 1: for p > | takes the form /? > 1, which always holds, and for 
/7 < I it reads p > 3 — j. These two conditions are: p > 2 - |/?. D 

Please recall, that by point 2. merely from existence, i.e for p - p, one has w < /? - |. 

3. Conclusion. 

3.1. Note on the rigorous estimates. 

The above computations are formal. To perform them rigorously, transform the considered prob- 
lem using Steklov averages with respect to time and use finite differences instead of differentiat- 
ing it with respect to space. This procedure has been presented in DiBenedetto |i^, DiBenedetto 
and Friedman 191] for homogenous systems. In our case we need to deal additionally with a 
quasilinear right-hand-side, for which the testing function u\ |Vm|'?7^ may not be admissible. In 
order to begin iteration, we need to have: s^ + w + I < p and to perform it at the i-th step: 
Si + w + \ < Si + M. However the latter inequality holds from the definition of M, the former 
may prove sometimes troublesome. In such cases one can resort to testing with F„{u[ IVmI't/^), 
where F,j{x) is a Lipschitz truncation at the level n. As the estimates are valid for every « we 
can proceed as before. Observe however, that we do not encounter these difficulties during com- 
putations for Theorem[T] For additional rigorous treatment (especially for s nonpositive, consult 
Choe and Bae ||7[] and reference thereirQ as well as Choe lalsll). 



^Observe, however, that in Choe and Bae (3 there are allowed .so > -2, which seems to be incorrect as far as 
i £ (-2,-1] are concerned. 
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3.2. Further research. 

The most obvious possible generalizations to the result is to allow for faster growths of the right- 
hand-side at the cost of additional assumptions, especially as some of them, like boundedness of 
the solution, appear naturally in the existence theory. It would be interesting to obtain a general 
result for critically growing right-hand-side (i.e. like 1 -i-|Vm|'0, with some smallness assumption, 
which would generalize the classical results for the heat system, mentioned in the introduction. 
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